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In biochemical systems the Michaelis-Menten (MM) scheme is one of the best-known models of the enzyme- 
catalyzed kinetics. In the academic literature the MM approximation has been thoroughly studied in the context 
of differential equation models. At the level of the cell, however, molecular fluctuations have many important 
consequences, and thus, a stochastic investigation of the MM scheme is often necessary. In their work Barik et 
al. [Biophysical Journal, 95, 3563-3574, (2008)] presented a stochastic approximation of the MM scheme. They 
suggested a substitution of the propensity function in the reduced master equation with the total quasi-steady- 
state approximation (tQSSA) rate. The justification of the substitution, however, was provided for a special case 
only and did not cover the whole parameter domain of the tQSSA. In this manuscript we present a derivation of 
the stochastic tQSSA that is valid for the entire tQSSA parameter domain. 


I. INTRODUCTION 


concentration of the enzyme-substrate complex ( C ). TQSSA 
implies: 


Enzyme-catalized reactions are frequently encountered in bi¬ 
ology and are modeled with the MM scheme: 


E + S % CHE + P. (1) 

fc-i 

Here k\ and [ are reaction constants for the complex (C) 
formation and dissociation respectively, and k 2 is the catal¬ 
ysis rate constant. Rate equations for the three elementary 
reactions in the scheme are presented through the following 
ODEs: 


dS_ 

dt 

C(S) 


~k 2 C{S ), (6) 

2E 0 S 

E 0 + K m + S+ s/{E 0 + K m + S ) 2 - 4 E 0 & 

(7) 


The uniform validity domain of the tQSSA is k-\ k 2 . 

However, any one of the following domains is also sufficient 
for the validiy of the tQSSA 0 : 


— — — kiS(E 0 — C) + fc-iC, (2) 

dt 

dC 

= —(fc_ 1 + k 2 )C + kiS(E 0 — C). (3) 

dt 


E 0 < Km + So, (8) 

S 0 -C E 0 + K m and fc_i > k 2 , (9) 

E 0 > S 0 and E 0 » I<m and k-\ -C k 2 . (10) 


Here Eq is the initial concentration of enzyme molecules, and 
the law of mass action is used to write enzyme E as (Eq — C ), 
assuming (Eq, So, 0,0) as the inital condition. 

In 0 ] Segel derives the deterministic quasi-steady-state ap¬ 
proximation (QSSA) with the validity condition 

fio«So + A M , (4) 

where Km = (fc-i + k 2 )/ki is the MM constant. The ap¬ 
proximate rate law reads: 


In order to investigate the stochastic behavior of the MM 
scheme, a chemical master equation is introduced. This leads 
to a reduced master equation for substrate or total substrate 
molecules, providing an approximation to the exact solu¬ 
tion. 

In Jit] Barik et al. introduce a general method of dealing with 
the stochastic tQSSA and derive the following reduced master 
equation for the MM scheme: 


dS __ k 2 E 0 S 

dt ~ K m + S’ ; 

Tzafriri 0 extends the results obtained by Segel 0 by sug¬ 
gesting the method called the total QSSA (tQSSA). It ap¬ 
plies to reactions where the concentration of the total sub¬ 
strate (S = S + C) changes on a slower timescale than the 
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^ = k 2 (c\s + 1 )P{s + 1; t) - k 2 (c\s)P(s; t ), 
at 

(id 

where (c|s) is the conditional expectation of the number of 
complex molecules at the quasi-steady-state. The authors then 
substitute (c|s) with the quasi-steady-state value of the com¬ 
plex c(s) from the tQSSA obtained by Tzafriri 0 . To support 
the validity of the method, they derive a recurrence relation for 
the quasi-steady-state probability distribution P(c|s): 
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P(c|s) 


f[ -jp— (^ ~ ^ + l)( e o —i + 1 ) p (Q| g )) 

i=i fc-1 * 

0 < c < min(s, eo) (12) 


and compare (c|s) obtained from Equation (O with the 
tQSSA result for a specific choice of parameters. Seeing an 
almost exact matching, the authors conclude that (c|s) ac¬ 
quired from Equation (fTIli ubiquitously substitutes the tQSSA 
result given by Equation (|7}. However, the mere fact that 
Equation (ITH i does not contain the Michaelis-Menten constant 
Km while Equation (|7} does, already suggests an inconsis¬ 
tency. 

In what follows we provide a derivation of the stochastic 
tQSSA for the MM scheme, based on the underlying assump¬ 
tions of the tQSSA. Unlike the argument presented by Barik et 
al. 1 which holds for a limited choice of parameters, our ap¬ 
proach holds for the entire tQSSA parameter domain. 


II. STOCHASTIC TQSSA 



FIG. 1: A comparison of the dynamics of the substrate molecules 
obtained from the Gillespie simulation of the scheme 0 with that 
obtained from the stochastic tQSSA given by Equation 1221 . The 
solid line represents the exact mean and the dashed lines represent 
one standard deviation away from the mean. The starred lines rep¬ 
resent the tQSSA mean and ±1 standard deviation away from the 
mean. Parameters used: k\ = 0.1s -1 , k- 1 = 5s -1 , k 2 = 5s -1 , 
.Eo=10, So= 10. TQSSA condition ((Sjl is met: Eq So + Km 


We begin by introducing the chemical master equation for the Multiplying both sides of Equation £[6]) by c and summing 
total substrate s = s + c and the complex c. over c, we obtain: 


= ~ [ki(s-c)(e 0 -c) + (k - 1 + k 2 )c]P(s,c-,t) 

+ ki(s — c + l)(e 0 — c + l)P(s, c — 1; t) 

+ k-i(c + l)P(s, c + 1; f) 

+ k 2 (c + l)P(s + 1, c + 1; t). (13) 

Summing Equation (IT3l > over c, we obtain: 


d(c\s\ t) 
dt 


where 


ki (se 0 - (e 0 + s)(c|s; t) + (c|s; t) 2 + a^ s (t)) 
(k -i + k 2 )(c\s;t), 

(17) 


dP(s; t ) 
dt 


where 


-k 2 (c\s; t)P(s ; t) + k 2 (c\s + 1; t)P(s + 1; t), 

(14) 


(*) = - ( c l s; *» 2p ( c l s ; *)■ ( 18 ) 

c 


Applying the tQSSA, we set the right side of Equation (IT71) 
equal to zero, obtaining: 


(c|s;f) = ^2cP(c\s;t). 

C 


(15) 


(c|s; t) 2 - r(s)(c\s;t) + e 0 s( 1 + 




eo s 


= 0, (19) 


The evolution of the complex c for the given number of the 
total substrate s is governed by: 


where 


^^dt ^ = ~ [^i( s - c)(e 0 - c) + (fc_i + fc 2 )c]P(c|s;f) 

+ ki(s — c + l)(eo — c + l)P(c — l|s; t) 

+ (fc_i + k 2 )(c -(- l)P(c + l|s; t). (16) 


r(s) = K m + e 0 + s. (20) 

Typically cr 2 ,-(t) is on the same order of (c|s;£), which, in 
turn, does not exceed the initial number of enzyme molecules 
eo- Thus, considering also the total substrate s term in the 
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FIG. 2: The ratio a'f.^/eo's appearing in Equation d 1 as a function 
of s at three different times, obtained from the exact solution of the 
original system G3- Parameters used: k\ = 0.1s r , k -1 = 5s 1 , 
k 2 = 5s _1 , Eo=lO, So=10. TQSSA condition {8} is met: Eq -C 
So + K m 

denominator, we can safely neglect the ratio tr^-(f)/eoS in 
Equation ( IT9l ). (Fig. [2] depicts the ratio as a function of the 
total substrate s for three different times at the quasi-steady- 
state. As we can see, in all three cases it is much smaller than 
1 for all values of s .) Solving for (c|s; t) at the quasi-steady- 
state, we obtain: 


r(s) — y/r*J7)—^4eos 

(c| S ) « - 2 - 

2eos 

eo + Km + s + \/ (eo + Km + s) 2 — 4eos 

( 21 ) 

This expression matches exactly with the tQSSA result 0 
obtained by Tzafriri (2], meaning that the substitution of the 
propensity function with the tQSSA rate function in Equation 
(fl4l> is justified under the tQSSA conditions. This results in 
the reduced master equation for the total subtrate s: 

dP(S' 

—^— = ~k 2 {c\s)P(s-, t ) + k 2 (c\s + 1 )P(s + 1; t), 

( 22 ) 

where (c|s) depends only on the total substrate s and is given 
by Equation (|2IT >. Fig. Q]compares the dynamics of the num¬ 
ber and of the variance of substrate modecules obtained from 
the Gillespie simulation with those obtained from the stochas¬ 
tic tQSSA given by Equation (f22l> . 

Applying the tQSSA (P(c\s;t) « 0) on Equation ( fl6l ). we 
can also derive a recurrence relation for the quasi-steady-state 
probability distribution P(c|s): 


FIG. 3: A comparison of the conditional expectation (c|s) obtained 
from the recursive relation and Equation S3 with that obtained 
from the exact solution of the original system G3 at three different 
quasi-steady-state moments: t = Is, 3s, 7s. The dashed line rep¬ 
resents Equation d, the solid line - Equation CD- and the cir¬ 
cles lines - the exact solution. Parameters used: k\ = 0.1s -1 , 
k-x = 5s -1 , fe 2 = 5s -1 , E 0 =W, 50=10. TQSSA condition ® 
is met: E 0 < So + Km 


j=(c|») = ft 2- (5 - i + 1) ( e °-‘ + 1) p ( o|»), 

. "t A m % 

2=1 

0 < c < min(s, eo). (23) 

Notice that the relation d23l) is very similar to the relation 
(fl2l i. obtained by Barik et al. {H]. The only difference is in 
the front coefficient, which in our case is 1 /ATm> whereas in 
Equation (fl2l) it is fci/fc_i. That is why the proof presented 
by Barik et al. is valid only when Km ~ fc_i//ci, equiva¬ 
lent to k 2 /k-i « 1. Fig. [3] demonstrates that (c|s) given 
by Equation (fTil) matches with (c|s; t) obtained from the ex¬ 
act solution at three different quasi-steady-state moments (t 
= Is, 3s, 7s) much more closely than (c|s) obtained by the 
recursive relation (IT2l > suggested by Barik et al. This means 
that the method of Barik et al. is indeed limited to the case 
k-x » fc 2 - 


in. DISCUSSION 

In this paper we presented a derivation of the stochastic 
tQSSA for the MM scheme that was valid for the entire 
parameter range of the tQSSA. Substrate dynamics obtained 
from the reduced master equation was in perfect agreement 
with the Gillespie simulation of the MM scheme. We also 
showed that the method suggested by Barik et al. was 
limited to the special case k-\ <C k 2 . To our knowledge, our 
derivation of the stochastic tQSSA is a novel contribution. 














All code used in this manuscript is available as open source on approximations to the stochastic MM scheme. 

https://github.com/galstyan/stochasticTQSSA, 
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We would also like to mention that an alternative approach of 
reducing the stochastic descriptions of biochemical networks 
was recently suggested by Thomas et al. 0. They used the 
method called slow-scale linear noise approximation (ssLNA) 
to evaluate the intrinstic noise present in enzyme networks. In 
the limit of large system size, the ssLNA appears to give better 
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